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Non-linear electrodynamics of a ring-shaped Andreev inter- 
ferometer (superconductor-normal conductor-superconductor 
hybrid structure) inductively coupled to a circuit of the dissi- 
pative current is investigated. The current-voltage character- 
istics (CVC) is demonstrated to be a series of loops with sev- 
eral branches intersecting in the CVC origin. The sensitivity 
of the transport current to a change of the applied external 
magnetic flux can be comparable to the one of the conven- 
tional SQUID's. Spontaneous arising of coupled non-linear 
oscillations of the transport current, the Josephson current 
and the magnetic flux in Andreev interferometers are also pre- 
dicted and investigated. The frequency of these oscillations 
can be varied in a wide range, while the maximal frequency 
can reach LO ma x ~ 10 12 sec -1 . 

Recent years much attention has been paid to the 
charge transport in mesoscopic systems which combine 
normal conductor (N) and superconductor (S) elements 
(for review papers see, e.g., [2,3] and references there). 
In such a hybrid structure superconducting correlations 
penetrate into the normal conductor within mesoscopic 
distances changing its transport properties there. This 
quantum effect is most pronounced in S-N-S structures 
(" Andreev interferometers" ) in which the quantum in- 
terference gives rise to a high sensitivity of these super- 
conducting correlations to the phase difference between 
the superconductors <p. It results, in particular, in os- 
cillations of the Josephson current J s — J s (<p) with a 
change of <p, the amplitude of the oscillations [3] being 
J c ~ N±evp/L]y (N± = Sn/X%, \f — h/PF - the Fermi 
wave length, pp - the electron Fermi momentum, e - thc 
electron charge, vf - the Fermi velocity, Sn and Ljv are 
the cross-section area and the length of the normal con- 
ductor, respectively) . 

If the Andreev interferometer is coupled to reservoirs 
of normal electrons, the conductance G(tp) and hence the 
dissipative current, Jd = G(ip)V, between the reservoirs 
(V is the bias voltage) are also phase-dependent, the am- 
plitude of the conductance oscillations being proportional 
to N±e 2 /Ti (in many experiments N± ~ 10 2 -f- 10 5 ). 

In a ring-shaped geometry of the S-N-S structure, the 
phase difference between the superconductors is con- 
trolled by the magnetic flux $ = HS r (H is the mag- 
netic field thrcading the ring, S r is the area of the ring): 
<p = — 27r<í>/$o where the flux quantum <í>o = nhc/e, c 
-the light velocity. This allows to control the Joseph- 
son and the dissipative currents with a change of the 
extcrnal magnetic field H ext that is favorable for various 
applications. On the other hand, the Josephson current 



(which depends on the magnetic flux in the ring) cre- 
ates its proper magnetic field which modifies this flux. 
As a result, for a ring of a large enough self-inductance, 
the relation between the magnetic flux $ and the ex- 
ternal magnetic field H ext turned out to be highly non- 
linear with hysteresis loops in the dependence of $ on 
H ext as shown in Fig.2 (see, e.g., [4]). This self-inductive 
property of the Josephson current has been used for cre- 
ation of Supcrconducting Quantum Interference Devices 
(SQUID) which have found extensive applications for ex- 
tremely precise measurements [4] . 




FIG. 1. Superconductor-normal conductor-superconductor 
structure of the Andreev interferometer type to which a volt- 
age drop V is applied. Thick lines indicate potential barriers 
at normal conductor - superconductor interfaces 1 and 2, and 
between the lead and the normal section of the interferometer. 

A dissipative current fiow through the S-N-S ring cre- 
ates a new situation in which the Josephson and dissi- 
pative currents are coupled together through their joint 
influence on the magnetic flux inside the S-N-S ring that, 
in its turn, affects the both currents themselves. 

The aim of this paper is to show that the inductive 
interaction between the dissipative and Josephson cur- 
rents results in a complicated loop-shaped form of thc 
current-voltage characteristics (CVC) of S-N-S structures 
schematically presented in Fig.l. In the general case sev- 
eral branches of the CVC (which correspond to different 
values of the magnetic flux <I> inside the ring) intersect 
in the origin of CVC ( J¿l = 0, V = 0) as is shown in the 
insertion of Fig.(4). Stable non-linear periodic in time 
oscillations of the dissipative current J¿, the Josephson 
current J s and the magnetic flux <I> inside the ring (and 
hence the phase diffcrence <p) are predicted and investi- 



f 



gated. The form of the loop-shapcd CVC and the fre- 
quency of thesc oscillations are shown to be extremely 
sensitive to the applied external magnetic field H ext and 
the bias voltage V. 

The system under consideration is a ring-shaped supcr- 
conductor which has a section of a normal conductor (see 
Fig.l). A steady bias voltage V is applied between the 
normal and the superconducting sections of the ring that 
provides a dissipative current Jd in this circuit. There is 
a potential barrier of a low transparency t r between the 
normal conductor and the lead. The electron transport 
may be in the ballistic or in the diffusive regimes. 
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FIG. 2. Dependence of the magnetic flux threading the 
ring, <í>, on the applied flux & ext = H ext S r . 

1. Loop-shaped current-voltage characteristics. The 
total flux of the magnetic field threading the ring, •!>, is 
created by the external magnetic field H ext , the propcr 
magnetic fields of the Josephson current J s and the dis- 
sipative current J<¿, and hence can be written as 



<&ext + —Js 

C 



C 



12 



Jd 



(1) 



where <& e xt — H ext S r is the magnetic flux of the exter- 
nal magnetic field H ext , C r is the self-inductance of the 
ring, C\2 is the mutual inductance of the ring and the 
dissipative current circuit. 

Eq.(l) togcther with the conventional relation Jd = 
G(ip)V define a parametric form of the CVC: 



Jd = Jd(<P) = -jr- [~~^(<Pext 



V 



Mf) 
G (<p) : 



Pext 



ip) - C r J s (ip) 

27T$ext/$0 



(2) 



One of the distinguishing features of the CVC of the 
system under consideration is its "many-fold degener- 
acy" , that is several branches of the CVC (which corre- 
spond to different values of the magnetic flux threading 



the ring <3>) can intersect in its origin (V = 0, J<¡ = 0). 
As is seen from Eq.(2) at J¿ = and Fig.2, the "degener- 
acy factor" is equal to the number of intersections of the 
vcrtical line & ext = const and the curve <3> = <3>(<3> ea;t ). 

Other key features of the CVC can be found if one con- 
sider the differential resistance dV/dJd which is readily 
obtained from Eq.(2): 



dV 
dJ d 



G(<p) 



1 



2tt C 12 J d G'(p) \ 

c$ G(cp) A(<p) ) v=7?(Jd) 



(3) 



Here and below the sign "prime" means the derivative 
with respect to <p; the function A(<p) [5] is 



Mv) = i 



2nC r 
c<3>o 



(4) 



All quantities in Eq.(3) are taken at (p = Tp(Jd) which is 
a solution of the first equation in Eq.(2). Using Eq.(2) 
one finds 



dTp/dJd 



2tt£i 



1 



c<3>o A(ip) 



(5) 



If C r > C er ] = c<P /(2Trmax{\J' s \}) there are values 
of <p at which A((p) = 0. Thcrcforc, in this foi- 
lows from Eq.(3) the CVC Jd(V) inevitably has points at 
which thc diffcrential resistance dV/dJd goes to infinity 
changing its sign there. On the other hand, the deriva- 
tive G' also changes its sign with a changc of ip that can 
provide points at which dV/dJd goes to zero changing 
its sign at thcm. The consequtive order of this pecu- 
liarities with an increase of Jd (that is thc form of the 
CVC) depends on the relative positions of the maxima 
and minima of G(ip) and J s (p). 

For the case of a low transparency of the potential 
barrier i r Cl the conductance G((p) has maxima at odd 
numbcrs of n [2] while thc Joscphson currcnt J s = 
at them and has a maximum and a minimum in their 
vicinities [3] as schematically shown in Fig.3. This gen- 
eral information together with Eq.(3 - 5) would suffice to 
find the CVC to be loop-shaped if C r > C er ^ . In order to 
see it let us start from the external magnetic field which 
corresponds to Tp = and hence G' — 0, J' > and 
therefore dV/dJd > according to Eq.(3). As at this 
point dp/dJd < the phase difference Tp decreases with 
an increase of current Jd (see Eq.(5)). Therefore, with 
an increase of the current G' becomes negative and after 
passing the minimum of J s (<p), AÍTp) starts to decrease 
because J' < now (see Fig.3), and when A is small 
enough (but positive yet), the differential resistance be- 
comes dV/dJd = 0. With a further increase of the cur- 
rent AÍTp) — > +0 while G' is yet negative, and hence 
dV/dJd — > — oo. Whcn A(Tp) = —0 thc differential re- 
sistance dV/dJd = +oo. In order to follow this second 
branch of CVC one should decrease the current because 
dTp/dJd > (scc Eq.(5)). Pursuing such a reasoning 
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one easily finds the current-voltage characteristic to be 
a serics of loops which touch the lines J¿ — G m i n V and 
■h = GmaxV (in our case G min = G(0), G max = G(ir)); 
the number of loops intersecting in the origin of the CVC 
increases with an increase of the ring inductance. Such 
a loop-shaped CVC takes place for the both cases of a 
diffusive or a ballistic normal conductor of the S-N-S ring- 
shaped interferometer. 
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FIG. 3. Typical dependences of the Josephson current and 
the conductance on the superconductor phase difference <p 
for the case of a low transparency (t r -C 1) of the potential 
barrier between the normal section of the sample and the lead 
(J c = N ± ev F /L N and G = N ± e 2 /{irh)). 

An example of such a current-voltage characteristics 
is presented in Fig.4 for the case of a ballistic normal 
section of thc Andreev interferometer in the presence of 
potential barriers at the N/S interfaces. In this case the 
dependences of the conductance G(<p) [7] and the Joseph- 
son current J s (<p) can be written as follows : 



G(<p) = 



t 2 Go 



(l + |r^^ 2) |cos^ + í2/2) -|r<¿ 



(6) 
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J s (<p) = N ± -^\ryry\smip x 
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N 



Ja 



g 2irsm(p-((f>) 



arctan 



' sm<p-(<fi) 



1 — e -2 *,- cos <P-(4>) 



Here G = N ± e 2 /(wh), 



N ± 



Vñ 



(7) 



VF 



and Vf¡ is the electron velocity in the n-th transverse 



mode; 



„(1.2) 
' N 



are the normal reflection amplitudes at 



N/S interfaces 1 and 2 whilc 

[8]); cos(^_(0) = 
sin <p-(4>) = y/í 



cos< 



-r 7jv 

|Jl) r (2) 
\ r A r A 



= 1, (see 
cosí/? and 



COS 2 Íy9- 



((/)). Expanding Eq.(7) in 



(1,2) 



onc 



exp{(p_} and taking the limit — > 0, ¿ 
reduces it to the wcll-known expression for the Joscph- 
son current in a 3D SNS junction [9]. 

Inserting Eq.(6,7) in Eq.(2) and performing numerical 



calculations one obtains the current-voltage characteris- 
tics shown in Fig.4. 
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FIG. 4. Current - voltage characteristics of an SNS ring 
for the external flux $ elt /<3?o = (21 + 1)tt, / = 0, ±1, ±2, 



t T = 0.1, = 0.2, t N ' = 0.25, and C r = 2.5c$ /J c ; 
J = (C r / '£12) J c , Vo = 2hv F /(ireL N ). The insertion shows a 
zoomed vicinity of the CVC origin. 

2. Electromagnetic self-oscillations. A change of the 
magnetic flux thrcading a supcrconductor-normal con- 
ductor ring results in the following current in the ring: 

1 d^> 
aR~dt 

where R is the resistance of the normal section of the 
S-N-S ring [10]. Using <p = -2tt$/$ onc sces Eq.(8) to 
be the equation of the ac Josephson effect. 

Taking into account the self-inductances of the ring C r 
and thc transport current circuit C\\ togcther with their 
mutual inductances £12, after using Eq.(l) (in which J s 
should be changed to J r ) and Eq.(8) one gets a set of 
equations that describes the time evolution of the trans- 
port current and the superconductor phase difference as 
follows: 



,(2) 



J r = J s (<p) 



(8) 



C 12 J' s (<p) d<p 



+ 



Jd 



C^\dJd_ C\ 2 <&Q d 2 <p 
c 2 dt 2irc 3 R dt 2 
$ £ r dw c$ , . 

r=— + ——<p + C r J s (<p) + £12 Jd = -C^ext 

2-KcR dt Z7T 



dt G(<p) 



= V: 



(9) 



The static solutions of Eq.(9) give the CVC (Eq.(2)). 

In ordcr to investigate time-evolution of the system it 
is convenient to eliminate J<¿ from the set of equations 
Eq.(9) and get the following closed equation for <p(t)\ 



(10) 



where C e ff = (C r C\\ — C 2 2 )/C\\). Eq.(10) is the equa- 
tion for a non-linear oscillator under a "friction" 
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70) = 
and a "force 



Cu RG(<p) 



(11) 



<P - <Pext + 



CnG(tp) 
2irC r 



J s (<p) 



2-kC x 
c$ 



(12) 



For low values of Cn the "friction" 7 > at any valuc 
of <¿j and hence, starting from any initial state, the system 
approaches one of its static states determined by Eq.(l) 
(in which Jd = G(<p)V). With an increase of £n, the 
"friction" j(<p) becomes negative in a certain range of <p 
and the static state can become unstable. Investigations 
of the stability of the static solutions of Eq.(lO) <p = <p s t 
show that the critical value of £n is determined by the 
condition j(<p s t) = 0, that is 



B 



BRG(<p st ) : 
2nC eff 



1 + 



c$ 



JÍM > 



(13) 



The Poincare method [11] shows that in the plane (<p, <p) 
a stable limit cycle arises if 



£11 > C cr ; 6 = 1 



-JÁPst) + —z—G'fait) > 



c$ 



c$ 



These inequalities (together with thc onc in Eq.(13)) can 
be satisficd if only J' s (<p st ) < and G'(<p st ) > 0. From 
here and Eq.(3,4) it follows that the stable limit cycle 
(that is non-linear periodic time-oscillations of <p(t) and 
<p(t)) can only be on those branches of thc CVC with 
the negative diffcrential resistance dV/DJd < that are 
close to the Jd-axis in Fig.2 (and on thc loop shown in 
the insertion of the figure). If < £11 — £ cr -C C cr the 
frequency of these oscillations is 



u>o = 



bc 4 R 



(CuC r - C\ 2 ) G(<p st 



This frequency can be variated in a wide range, and esti- 
mations show the maximal frcqucncy u> max ~ 10 12 sec -1 . 

In conclusion, for both the diffusive and ballistic cases 
I have shown that the current-voltage characteristics of 
the Andreev ring-shapcd intcrfcrometer Jd(V) is a series 
of loops with several branches intersecting in its origin, 
and hence it has sections with a negative differential re- 
sistance dV/dJd if the ring inductance is large enough. 
These properties of the CVC are robust as the negative 
differential resistance and the intersection of the several 
branchcs in the CVC origin appear at any value of the 
mutual inductance, £12, between the ring and the dis- 
sipative transport circuit (which can be as small as is 



wished) as soon as the self-inductance of the ring C r is 
large enough to provide A < (see Eq.(3,4); of course, 
such a CVC can not be observed experimentally if the 
mutual inductance C12 is unreasonably small. This in- 
cquality (A < 0) is the same as the one needed for the 
conventional SQUID functioning. What's is more, ma- 
nipulations with the ring self-inductancc C r and the ex- 
ternal magnetic fiux $ ext allow to get CVC loops (and 
therefore the hysteresis loops of the CVC) as small as 
wanted, especially it concerns the smallest loop around 
the CVC origin (see Insertion in Fig.4), and one can ob- 
tain the sensitivity of the transport current Jd (or the 
applied voltage V) to a change of the external magnetic 
flux comparable to the conventional SQUID sensitivity. 

Spontaneous arising of coupled non-linear oscillations 
of the transport currcnt Jd(t), the Josephson current 
J s (t) and the flux <E>(í) in Andreev interferometers have 
been also predicted and investigated. The frequency of 
the oscillations u> can be varied in a wide range, and thc 
maximal frcquency can reaches uj max ~ 10 12 sec -1 . 
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